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ABSTRACT 


An  improved  model  for  a  once-through  counter-cross-flow 
waste  heat  recovery  unit  with  a  segmented  fin-tube  arrange¬ 
ment  was  developed.  The  model  was  coded  in  FORTRAII  IV 
computer  language  for  use  on  an  IBM-3033  computer  system/ 
and  was  tested  for  various  conditions  of  uniform  and  non- 
uniform  gas  glow  distributions.  Additional  parameters 
varied  were  heat  exchanger  size,  stecun  outlet  pressure, 
steam  outlet  temperatxire  and  gas  flow  rate.  The  effect  of 
each  variation  on  WHRU  performance  was  evaluated  and  the 
interrelationship  discussed.  Results  indicate  the  ideal 
design  is  a  high  steam  temperature,  low  operating  pressure, 
flow  distribution  controlled  heat  exchanger. 
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p  -  Density  (Ibm/ft) 

-  Density  of  Saturated  Water  (Ibm/ft) 

-  Density  of  Saturated  Vapor  (Ibm/ft) 


I .  INTRODUCTION 


A.  BACKGROUND 

Since  the  introduction  of  the  DD-963  class  ships  into 
the  fleet,  the  U.S.  Navy  has  become  firmly  committed  to  gas 
turbine  propulsion  plants.  Currently,  besides  the  DD-963 
class,  there  are  four  additional  ship  classes  which  are,  or 
will  be,  gas  turbine  powered.  These  are  the  Patrol  Hydro¬ 
foil  Missile  (PHM's) ,  the  Oliver  Hazard  Perry  FFG-7  class, 
the  Ticonderoga  CG-47  class  and  the  soon  to  be  built  DDG-51 
class. 

This  commitment  to  gas  turbine  propulsion  has  stimulated 
great  interest  in  the  concept  of  converting  the  waste  heat 
from  the  gas  turbine  engines  into  useful  shaft  power.  In 
December  1979  the  Naval  Sea  Systems  Command  awarded  contracts 
for  development  of  a  conceptual  design  of  a  new  generation 
of  waste  heat  recovery  systems.  The  goal  of  the  design  was 
to  provide  a  system  which  would  realize  a  25  percent  minimum 
improvement  in  ship  propulsion  fuel  consumption  and  increase 
the  ship's  operating  range  by  1000  nautical  miles.  The  out¬ 
growth  of  the  effort  has  been  the  RAnkine  Cycle  Energy  Recovery 
(RACER)  system.  As  conceived,  the  RACER  system  will  be  an 
unfired  waste  heat  recovery  system  designed  to  convert  waste 
heat  from  the  exhaust  of  the  main  propulsion  gas  turbine 
engines  into  useful  shaft  horsepower.  Such  a  system  has 

15 


coimnonly  been  called  a  COGAS  or  combined  Gas  And  Steam  system. 
The  hot  exhaust  gases  of  the  main  propulsion  turbine  are 
passed  through  a  Waste  Heat  Recovery  Unit  (WHRU)  ,  where  steeun 
of  the  desired  temperature  and  pressure  is  generated.  The 
steam  is  used  to  drive  a  steam  turbine  which  operates  in 
parallel  with  the  propulsion  gas  turbine  through  a  common 
reduction  gear. 

Numerous  papers  have  been  written  concerning  the  RACER 
system.  Halkola/  Cambell  and  Jung  [Ref.  1]  discussed  the 
conceptual  design  for  the  system.  Mattson  [Ref.  2]  presented 
a  plan  for  designing  reliability  and  maintainability  into 
the  RACER  system  and  Combs  [Ref.  3]  developed  a  computer  model 
of  a  COGAS  system  which  designs  a  waste  heat  recovery  unit 
and  allows  comparison  of  system  performance  by  variation  of 
system  parameters. 

References  1,  2  and  3  all  agree  that  the  desired  goal  of 
fuel  conservation  must  not  come  as  a  result  of  increased 
maintenance,  mempower,  or  main  propulsion  system  down-time. 

The  design  philosophy  dictates  that  the  system  have  high 
reliability,  be  inherently  safe,  have  low  maintenance  require¬ 
ments  and  be  essentially  automatic  in  operation. 

Reference  1  states  that  the  waste  heat  recovery  unit 
paces  the  design  of  the  system.  The  majority  of  effort  by 
Combs  [Ref.  3]  was  to  create  a  WHRU  design  model  which  would 
simulate  the  performance  characteristics  of  a  physical  boiler. 
Since  great  expense  is  involved  in  the  design,  fabrication 
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£Uid  operation  of  prototype  systems,  it  is  of  primary  impor¬ 
tance  to  achieve  a  realistic  mathematical  simulation  first. 

To  this  end.  Combs  devised  a  model  which  designs  a  WHRU 
when  given  feedwater  inlet  temperature,  steam  outlet  tempera¬ 
ture,  operating  pressure,  inlet  and  outlet  gas  conditions 
and  gas  flow  rate. 

In  this  study,  the  initial  model  of  Combs  was  considered 
and  an  improved  computer  model  of  the  waste  heat  recovery 
unit  was  developed.  A  once-through  counter-cross-flow  heat 
exchanger  was  selected  because  of  the  advantages  of  lower 
initial  cost,  faster  response,  compactness  and  low  operating 
cost.  The  improvement  in  the  model  is  that  the  heat  exchanger 
is  discretized  into  individual  "nodes”  much  in  the  manner 
done  by  Shu  emd  Kuo  [Ref.  4].  This  procedure  treats  smaller 
physical  pieces  of  the  overall  heat  exchanger  as  individual 
small  heat  exchangers  to  attain  a  better  understanding  of 
the  temperature  distributions  within  the  WHRU.  The  improved 
model  allows  an  analysis  of  the  effect  of  non-uniform  velocity 
distributions  on  the  performance  of  the  WHRU,  and  ultimately, 
the  RACER  system. 

B.  OBJECTIVES 

There  were  five  main  objectives  for  this  thesis. 

1.  Devise  an  improved  computer  model  for  the  once-through 
WHRU.  This  model  was  to  predict  the  performance  of  various 
size  WHRU's  given  feedwater  inlet  temperature,  steam  outlet 


temperature,  WHRU  operating  pressure,  inlet  gas  mass  flow 
distribution  and  inlet  gas  temperature  distribution. 

2.  Enable  the  designer  to  estcUolish  a  design  WHRU  size 
based  on  the  performance  of  the  model  with  various  input 
conditions . 

3.  Evaluate  the  design  WHRU  utilizing  various  uniform 
and  non-uniform  inlet  gas  distributions  for  both  on-design 
and  off-design  conditions. 

4.  Re-evaluate  the  WHRU  utilizing  the  same  uniform  and 
non-uniform  gas  distributions  into  the  WHRU,  but  with 
different  steam  outlet  temperatures. 

5.  Provide  a  framework  for  future  studies  involving  a 
development  of  a  three  dimensional  model  of  a  WHRU  and  a 
detailed,  complete  COGAS  system  cycle  analysis  and  optimization. 


I I .  MODEL  DES CRI PT ION 


A.  OVERVIEW 

The  waste  heat  recovery  unit  chosen  was  a  once-through 
counter-cross-flow  heat  exchanger.  A  computer  program  was 
written  to  model  this  type  heat  exchanger  for  analysis  of 
the  effect  of  uniform  and  non-uniform  inlet  gas  mass  distri¬ 
butions.  Variation  of  other  parameters,  other  than  the  gas 
distribution,  provides  a  thorough  understanding  of  the  oper¬ 
ation  of  the  heat  exchanger  in  a  multitude  of  conditions. 

The  model  was  written  for  incorporation  into  a  RACER  system 
analysis,  but  it  is  not  limited  to  that  application. 

Figure  1  shows  a  three  dimensional  representation  of  a 
typical  layout  of  a  marine  gas  turbine  waste  heat  recovery 
unit.  The  unit  is  header ed  on  the  fluid  inlet  and  exit  and 
designed  to  be  installed  in  the  gas  turbine  exhaust  gas 
duct.  The  exhaust  gas  entering  the  bottom  of  the  WHRU  pro¬ 
vides  the  energy  for  transfer  to  the  fluid-side  of  the  system. 
The  fluid  (feedwater)  enters  at  the  upper  header  of  the  WHRU 
and  flows  alternately  back  and  forth  through  the  unit  until 
exiting  as  steam  at  the  lower  header.  The  U-bends  shown  in 
Figure  1  at  the  ends  of  the  tube  passes  allow  each  tube  to 
enter  and  exit  the  heat  exchanger  as  a  continuous  path  for 
fluid  flow.  As  a  tube  makes  one  horizontal  traverse  across 
the  WHRU  it  is  referred  to  as  a  pass.  Therefore,  the  number 


of  times  a  tiibe  passes  through  the  gas  flow  defines  the 
vertical  size  of  the  heat  exchanger  as  a  15,  20,  25,  etc., 
pass  heat  exchanger. 

A  gas  flow  distribution  is  depicted  as  entering  the  WHRU 
in  Figure  1.  The  gas  flow  is  shown  oriented  along  the  tube 
pass  as  a  non-uniform  distribution.  Additionally,  the  gas 
flow  may  be  non-uniform  in  the  direction  into  the  heat 
exchanger,  i.e.,  along  the  axis  of  the  headers  as  shown  in 
Figure  1.  This  two  dimensional  gas  flow  distribution  may  be 
non-uniform  in  both  mass  and  temperature.  The  nature  of  the 
heat  exchanger,  by  virtue  of  the  cross  flow  of  the  fluid, 
decrees  that,  even  if  the  mass  and  temperature  distributions 
entering  the  heat  exchanger  were  uniform  in  both  directions, 
these  distributions  can  become  distorted  and  non-uniform  as 
the  gas  passes  through  the  heat  exchanger  in  the  third 
dimension.  This  distortion  can  be  caused  by  the  variation 
of  the  heat  transfer  characteristics  along  the  tube.  The 
problem  of  analysis  becomes  more  complicated  if  cross  mixing 
of  the  gas  is  considered. 

The  model  of  Combs  [Ref.  3J  treated  the  gas  distribution 
as  a  uniform  and  constant  value  in  one  direction.  Combs 
assumed  the  mass  and  temperature  of  the  gas  along  the  axis 
of  the  tube  was  everywhere  the  same  for  a  complete  pass . 

The  model  yielded  results  which  were  close  to  the  results 
attained  by  other  studies  and  gave  a  reasonable  estimate  of 
the  overall  heat  transfer  characteristics  of  the  heat 
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exchanger.  The  model,  however,  could  not  predict  the  effect 
of  non-uniformity  of  the  gas  distribution  on  the  performance 
of  the  WHRU.  This  study  discretizes  the  dimension  along  the 
axis  of  the  tube  pass  and  enables  the  designer  to  evaluate 
the  effects  of  variations  in  mass  flow  rate  and  temperature. 
It  is  assumed,  however,  that  the  distributions  along  the 
axis  of  the  headers  are  all  the  same  as  the  outer  distribu¬ 
tion.  The  concept  of  discretization  and  "nodal"  arrangement 
is  shown  in  Figure  2.  The  representation  shows  a  matrix 
arrangement  where  the  number  of  rows  corresponds  to  the 
nvimber  of  passes  in  the  heat  exchanger,  and  the  number  of 
columns  corresponds  to  the  number  of  gas  paths  in  the 
discretization.  The  number  of  passes  is  a  user  defined  input 
to  the  model  but  the  number  of  gas  paths  is  fixed  in  the 
model  at  ten. 

Use  of  the  model  requires  that  the  user  provide  the 
following  conditions.  First,  the  number  of  passes  that  a 
particular  water/steam  tube  is  to  make  through  the  heat 
excheuiger  must  be  specified.  This  parameter  will  physically 
size  the  heat  exchanger  in  the  vertical  direction.  The  mass 
flow  rate  of  the  gas  turbine  exhaust  gas,  (m^) ,  for  each  of 
ten  exhaust  gas  paths  through  the  heat  exchanger  must  be 
entered.  For  a  uniform  distribution,  this  is  simply  the 
total  exhaust  gas  mass  flow  divided  by  ten.  A  simplifying 
assumption  of  the  model  is  that  there  is  no  cross  mixing  of 
flow  as  the  gas  passes  through  the  heat  exchanger;  therefore. 
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the  mass  entering  a  particular  gas  path  is  constant  through 
the  heat  exchanger. 

Next,  the  temperature  of  the  gas  entering  each  gas  path 
is  entered.  For  this  study  the  temperature  distribution  of 
the  gas  was  considered  uniform  for  all  runs;  however,  the 
model  is  not  constrained  to  uniform  temperature  distributions 

The  gas  was  assumed  to  be  at  atmospheric  pressure 
entering  the  heat  exchanger  and  gas  side  pressure  drops  were 
not  considered  in  the  calculations. 

The  user  next  enters  the  feedwater  inlet  temperature, 
the  steam  outlet  temperature  and  the  fluid-side  operating 
pressure.  These  parameters  are  used  in  the  initial  calcula¬ 
tions  to  determine  a  first  "guess”  mass  flow  rate  on  the 
fluid-side  of  the  heat  exchanger.  The  minimum  allowed 
exhaust  gas  temperature  is  entered.  This  value  is  a  user 
desired  constraint  and  is  specified  as  the  minimum  acceptable 
temperature  the  user  can  accept  to  prohibit  hot  side  corro¬ 
sion.  The  model  is  not  constrained  to  maintain  the  exhaust 
gas  temperature  above  the  user  specified  temperature. 

The  number  of  feet  in  length,  per  node,  is  entered. 

This  parauneter  is  based  on  the  physical  size  of  the  heat 
exchanger  being  analyzed  in  the  horizontal  direction.  For 
all  runs  in  this  study,  an  assumed  tube  row  length  of  ten 
feet  divided  by  ten  gas  paths  yielded  a  tube  length  per  node 
of  one  foot.  There  is  no  constraint  on  the  length  parameter, 
within  reason^d>le  limits. 


The  heat  exchanger  model  was  written  as  a  two  dimensional 
analysis,  i.e.,  all  parameters  in  the  third  dimension  were 
assumed  uniform.  The  next  parameter  entered  is  the  number 
of  tubes  per  row,  which  is  used  to  size  the  heat  exchanger 
in  the  third  dimension.  Figure  1  shows  a  three  dimensional 
representation  of  the  heat  exchanger  design.  Each  particular 
tube  enters  at  the  top  of  the  heat  exchanger  and  passes 
alternately  back  and  forth  below  itself  until  exiting  the 
heat  exchanger  at  the  bottom.  The  tubes  in  each  sxibseguent 
plane,  behind  the  first,  enter  and  exit  in  the  same  manner, 
but  alternate  slightly  up  and  down  to  give  a  staggered  tube 
arrangement.  The  nodes  are  numbered  in  standard  matrix  format 
with  the  upper  left  being  node  (1,1) ,  and  counting  across 
row  one  as  (1,2),  (1,3),  etc.,  to  (1,10);  then  node  (2,1) 
is  directly  below  node  (1,1) ,  and  so  on,  to  the  bottom  of 
the  matrix.  Figure  2  shows  the  nodal  arrangement  of  the 
heat  exchanger.  The  fluid  may  enter  from  either  side  of 
the  heat  exchanger.  If  the  heat  exchanger  has  an  even  number 
of  passes,  the  fluid  will  enter  on  the  same  side  as  the 
steam  exit,  cuid  on  the  opposite  side  for  an  odd  number  of 
passes.  Each  tube  is  continuous  from  entry  to  exit. 

The  final  user  entry  is  the  geometric  scaling  factor 
which  is  used  in  the  geometry  portion  of  the  model  to  scale 
all  the  physical  dimensions  of  the  tube,  with  the  exception 
of  the  length,  up  or  down. 


Two  additional  simplifications  were  assumed  in  the  model. 
First,  the  fluid  side  pressure  drop  was  not  considered. 
Lastly,  even  though  the  tube  passes  were  discretized  to 
determine  the  temperature  distributions,  the  tube  length 
within  a  node  still  has  a  significant  length  for  a  tempera¬ 
ture  gradient.  For  analysis,  the  gas  temperature  at  any 
node  was  assumed  to  be  uniform  across  the  length  of  the  node 
at  both  the  entry  and  exit. 


B.  INITIAL  CALCULATIONS 

The  WHRU  is  essentially  analyzed  in  three  sections; 
superheating,  boiling  and  heating.  Given  an  inlet  gas  distri 
bution,  water/steam  inlet  and  outlet  temperatures,  and  an 
operating  pressure,  there  is  a  unique  fluid-side  mass  flow 
rate  which  will  satisfy  all  boundary  conditions.  The  model 
attempts,  through  iteration,  to  determine  the  resultant 
fluid-side  mass  flow  rate.  The  analysis  is  also  greatly 
affected  by  the  number  of  passes  (i.e.,  the  physical  size) 
of  the  heat  exchanger.  There  are  clearly  combinations  of 
parameters  for  which  no  solution  is  possible,  or  the  solu¬ 
tion  is  trivial  in  that  the  mass  flow  rate  of  the  fluid  is 
so  small  as  to  be  meaningless. 


For  calculations,  variable  assignment  is  as  follows; 


(1)  mg(m,n) 


gas  flow  rate  through  gas  path  n  at  tube 
pass  m; 


(2)  T  (m,n) 
9 


gas  temperature  into  node  (m,n) .  The 
gas  temperature  out  of  node  (m,n)  will 
be  the  gas  temperature  into  node  (m-l,n) ; 


(3)  T^(ra,n)  =  fluid  temperature  into  node  (m,n) .  The 

fluid  temperature  out  node  (m,n)  will 
equal  the  fluid  temperature  into  node 
(m,n+l)  if  the  fluid  is  flowing  to  the 
right  (see  Fig.  3) ,  or  the  fluid 
temperature  out  of  node  (m,n)  will  equal 
the  fluid  inlet  temperature  of  node 
(m,n-l)  if  the  fluid  is  flowing  to  the 
left. 

At  the  end  of  a  fluid  pass  the  fluid 
out  of  a  node  will  equal  the  fluid 
temperature  into  the  node  directly  below, 
i.e.,  T^(m,n)  =  T^(m+l,n); 

• 

(4)  m^  =  fluid-side  mass  flow  rate 

Figures  3  and  4  are  representations  of  typical  nodes 
depicting  the  above  explanation  of  fluid  and  gas  flow. 

After  all  entries  are  made  the  pregram  initializes  by 
Slimming  the  individual  mg(NR,n)  elements  to  find  m^  total' 
then  the  Tg(NR,n)  elements  are  averaged.  The  T^(NR,n)  ele¬ 
ments  are  averaged  by  multiplying  each  element  by  the  corres¬ 
ponding  fraction  of  the  total  gas  mass  that  is  passing  through 
each  gas  path.  These  "weighted"  temperatures  are  summed  and 
the  result  is  divided  by  the  number  of  gas  paths.  The  result 
is  then  the  average  temperature  of  the  gas  into  the  heat 
exchanger.  An  average  gas  temperature  is  determined  using 
the  user  specified  minimum  gas  temperature  out  and  the  aver¬ 
age  gas  temperature  into  the  heat  exchanger.  Then; 


0  =  m  C  (T  -  T  ) 

“  g  total  pg  g  avg  in  g  avg  out 

where  C  is  evaluated  at  the  heat  exchanger  gas  bulk 

pg 

temperature.  Now,  a  first  "guess"  fluid  mass  flow  rate  is 


calculated  from  an  energy  balance 


where  h.  .  and  h- are  the  fluid  inlet  and  outlet 
£  in  £  out 

enthalpies  at  the  user  specified  feedwater  inlet  temperature 
and  steam  outlet  temperature. 

C .  GEOMETRY 

The  geometry  portion  of  the  model  is  nearly  identical  to 
that  used  by  Combs  [Ref.  3] .  A  fin-tube  with  helically 
curved  extended  surfaces  on  circular  tubes  is  assumed  for 
this  model.  The  fins  are  segmented  and  the  tubes  are 
configured  in  banks  of  one  row  each  with  the  rows  staggered. 
The  base  tube  surface  is  taken  from  Weirman,  Taborek  and 
Marner  [Ref.  5],  The  description  of  the  finned  tubes  is  as 
follows: 
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The  finned  tube  configuration  is  shown  in  Figure  5.  The 
tube  length,  L,  and  the  number  of  tubes  per  row,  are 

the  user  specified  entries.  The  tube  layout  is  shown  in 
Figure  6.  The  center  to  center  tube  spacing  in  the  trans¬ 
verse  direction  is  4.50  inches.  The  equilateral,  staggered 
tube  arrangement  used  in  this  model  leads  to  a  spacing 
normal  to  the  gas  flow,  s^^,  of  4.50  inches  and  a  spacing 
parallel  to  the  gas  flow,  s  ,  of  3.90  inches,  as  shown  in 

kr 

Figure  6.  The  heat  exchanger  height  is  established  by  the 

number  of  WHRU  passes  (rows)  and  the  tube  spacing  parallel 

to  the  gas  flow,  s  . 

P 

In  order  to  establish  the  minimum  gas  flow  cross-sectional 
area  the  total  "blocked"  frontal  area,  of  the  heat 
exchanger  must  be  calculated  from 

Ab  -  L  *  1.  Nj  J  tj  , 

and  the  minimum  gas  flow  area  is 


^rain  "  ~  \ 

where  is  the  total  frontal  area  of  the  heat  exchanger 
exposed  to  gas  flow.  The  total  inside  area  available  for 
heat  transfer  per  pass  is 


d.  L 


N 


t/r 
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The  gas-side  surface  area  available  per  pass  for  heat 
transfer  is  the  sum  of  the  fin  surface  area  and  the  bare 
tube  surface  area  per  pass.  The  fin  surface  area  per  tube 
is  calculated  from 


^fin 


N.LIN  (2  il  +  2  t,  £  +  t.)  +  Udl.  -  dj)] 

IS  cs  fc  sf  2  ro  o 


where  d^  =  d^  -  2£^.  The  bare  tube  area  per  tube  is 


*bt  =  »  <3o  ^  "  '*o  ’’f 


Therefore,  the  total  outside  area  per  pass  available  for 
heat  transfer  is 


^op  “  *^t/p^^fin  ■*'  \t^ 


The  cross  sectional  area  for  fluid  flow  is  calculated  from 


yd?  N  . 

4  1  t/p 


With  the  mass  flow  rates,  terminal  temperatures  and  heat 
exchanger  geometry  established,  the  remainder  of  the  model 
may  be  solved  by  nodal  analysis'  for  the  unique  fluid  mass 
flow  rate  and  the  location  of  the  beginning  of  the  boiling 
and  superheating  regimes. 


D.  SUPERHEATER  NODAL  ANALYSIS 

In  the  heat  exchanger  there  is  only  one  node  where  two 
conditions  are  initially  known,  and  that  is  at  the  node 
where  the  steam  is  exiting.  At  this  node,  both  the  gas  inlet 
mass  flow  rate  and  temperature  and  the  fluid  outlet  tempera¬ 
ture  have  been  specified  by  the  user.  Using  the  initial 
values  of  the  fluid  mass  flow  rate  from  the  estimated  overall 
energy  balance  and  the  desired  superheater  steam  outlet 
temperature,  the  heat  exchanger  is  analyzed  in  the  reverse 
direction  through  the  superheater,  through  the  boiler  and 
through  the  heater  to  arrive  at  the  feedwater  inlet.  For 
the  given  geometry  and  flow  conditions,  the  model  will  calcu¬ 
late  the  water  inlet  temperature  consistent  with  the  assumed 
estimated  mass  flow  rate.  If  the  temperature  arrived  at  by 
analysis  as  the  inlet  to  the  first  node  does  not  match  the 
user  specified  feedwater  inlet  temperature,  then  the  fluid- 
side  mass  flow  rate  is  adjusted  up  or  down  and  the  entire 
process  is  begun  again. 

To  begin,  at  the  final  node  the  gas  inlet  temperature 
and  mass  flow  rate  are  known.  Also,  the  fluid  outlet  tempera 
ture  and  a  first  "guess"  fluid  mass  flow  rate  are  known. 

It  is  then  assumed  that  the  bulk  (average)  fluid  and  gas 
temperatures  for  the  node  are  equal  to  these  known  tempera¬ 
tures.  All  fluid  and  gas  properties  for  a  node  are  evaluated 
at  the  bulk  temperatures. 

The  gas-side  Reynolds  number  for  the  node  is  calculated 
using  the  bulk  temperature.  With  the  Reynolds  number,  a 


j-£actor  is  obtained  from  a  polynomial  fit  to  the  data  for 
tube  layout  number  5  in  Ref.  5  (Fig.  6) .  By  definition, 
the  j-£actor  is  related  to  the  heat  transfer  coefficient, 
hg,  by 


j 


Pr 


2/3 


By  introducing 


St 


Nu 

Re  Pr  ' 


the  previous  expression  can  be  written  as 


j 


Nu 


9  9 


Pr 


2/3 


Nu 


and 


Nu  =  j  Re  Pry^ 
9  9 


where 


Nu 


h  d 


K 
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and  K  ®  the  thermal  conductivity  of  the  gas.  Therefore,  a 
9 

relationship  may  be  written  for  the  heat  transfer  coefficient 


as  follows: 


In  the  superheat  region,  the  water-side  heat  transfer 
coefficient  is  calculated  by  using  the  Dittus-Boelter 
correlation 


=  0.023  ^  Re°*®  Pr°‘'^ 


where  all  properties  are  obtained  at  the  bulk  temperature  of 
the  steam  in  the  node.  Using  the  tube  wall  resistance 
together  with  h^  cuid  h^,  the  overall  heat  transfer  coeffi¬ 
cient  for  the  node  Ceui  be  written  as 


U  . (m,n) 
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infdyap 


2itK  N 


w  t/p 


where : 


n,.  =  1  -  (1  -  T),) 


*fin 


f'  A 


op 


The  fin  efficiency,  is  calculated  from  the  expression 


tanh  ML 
ML 


where 


M  =  /h  P/KA 


and,  if  the  fin  is  approximated  by  a  set  of  rectangular 
strips  extending  from  the  tiibe  wall. 


L  = 


-  '*o 


Now,  the  cross-sectional  area.  A,  of  a  rectangular  strip 
may  be  written  as 


A  - 


and  the  perimeter  is 


P  =  2(Wg  +  t^)  . 


So,  using  known  geometric  parameters  and  the  thermal  conduc¬ 
tivity  of  the  fin  metal,  the  parameter  ML  may  be  expressed 
as 


where 


ML  =  C  /h“ 
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-  \\  P  > 
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Now, 


tanh  C/h~  A,. 
- 2_)  ^ 


C/h 


op 


[1  -  (1  - 


The  pass  effectiveness  is  calculated  using  the  effective- 
ness-NTU  method.  The  correlation  for  the  effectiveness  for 
a  single  pass,  cross  flow  heat  exchanger  with  both  fluids 
unmixed  is  taken  from  Incropera  and  DeWitt  [Ref.  6] 


e  =1-  exp[(;^)  (NTU)°*^^{exp[-C  (NTU)  °  -1}  ] 

p  C  r 


where : 
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C  .  is  the  lesser  of  C,  and  C  where 
min  f  g 


=£  =  V”* 


m„(m,n) 

g  pg  g 


and  C  is  the  greater  of  the  two.  All  fluid  and  gas 
max 

properties  are  again  evaluated  at  the  bulk  temperature. 

Two  equations  must  be  satisfied  to  find  the  solution  to 
a  node.  These  are  the  rate  equation 


Q  =  e  C  .  (T  .  -  T_  .  ) 

“  min  g  in  f  in 
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and  the  energy  equation 


Q  =  m,  C  -(T,  ^  . 

f  pf  f  out  f  in 


In  the  two  equations  the  only  unknown  is  Equating 

the  two  and  solving  for  yields 


,  _  min  g  in  f  pf  f  put 

min  f  pf 


Now,  Q  for  the  node  is 


Q  =  m-  C  ,(T-  ^  .  ) 

f  pf  f  out  f  in' 


and 
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Next,  new  bulk  temperatures  for  the  node  are  calculated  as 


"^g  bulk  ~  ^"^g  in  ^  "^g  out^ 


and 


"^f  bulk 


(T 


f  out 


+  T 


f  in 
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An  average  outside  tube  wall  temperature  is  calculated  from 


where  are  total  inside/outside  areas  for  the  node 

and  T  q/T are  the  new  node  bulk  temperatures .  The  expres 
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sion  for  T  _  is  derived  from  the  following  formulation  for 
wo 

heat  transfer  in  the  node 


which  reduces  to 
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The  gas-side  film  temperature  is 


This  gas-side  film  temperature  is  now  introduced  at  the 
beginning  of  the  calculations  for  the  node  as  the  new  gas 
bulk  temperature. 

The  assumed  bulk  temperatures,  at  the  beginning  of  the 
calculations,  were  obviously  too  high  since  the  gas  inlet 
and  fluid  outlet  temperatures  were  used.  This  assumption, 
however,  allowed  a  solution  for  the  node  to  be  found.  The 
model  stores  the  result  of  this  first  pass  and  returns  to 
the  beginning  utilizing  the  new  bulk  temperatures  to  calcu¬ 
late  a  new  solution  to  the  node.  At  the  end  of  the  second 
pass  through,  the  result  of  the  calculation  for  is 

compared  to  the  result  obtained  the  first  time.  If  T^ 
has  changed,  then  the  new  result  is  stored  and  the  process 
is  repeated.  When  no  longer  changes,  the  process  has 

converged  to  a  good  solution  for  the  node. 

A  check  is  made  to  determine  if  the  fluid  inlet  tempera¬ 
ture  is  less  than  or  equal  to  the  saturation  temperature  for 
the  pressure  specified.  If  the  temperature  is  higher  than 
the  saturation  temperature,  the  model  proceeds  to  the  next 
node  of  the  superheater  and  repeats  the  entire  process. 

If  the  fluid  temperature  equals  the  saturation  temperature 
superheating  has  theoretically  begun  at  the  interface  between 
this  and  the  final  node  of  the  boiling  portion  of  the  model. 

If  this  condition  exists,  the  model  proceeds  to  the  boiling 
nodal  analysis. 


As  mentioned  previously,  the  model  works  through  the 
heat  exchanger  in  the  direction  opposite  to  fluid  flow. 


When  it  is  determined  that  the  fluid  temperature  into  a  node 
is  below  the  saturation  temperature,  the  node  is  "split”  to 
determine  the  physical  location  in  the  node  where  super¬ 
heating  begins . 

To  split  the  node  the  m^  quantity  into  the  node  and  the 
physical  length,  L,  are  halved.  The  geometry  is  recomputed 
as  if  the  node  were  only  half  as  large.  The  solution 
process  is  reinitiated  and,  when  a  solution  is  achieved,  the 
saturation  temperature  check  is  redone.  If  the  temperature 

into  the  node  is  still  below  the  saturation  temperature,  m^ 

and  the  length  are  halved  again  and  the  process  is  repeated. 
This  procedure  is  continued  until  the  temperature  into  the 
node  is  found  to  equal  the  saturation  temperature  or  be 
above  it.  If  the  temperature  is  above  the  saturation  temp¬ 
erature,  then  the  ra^  and  L  parameters  are  changed  to  a  value 
halfway  between  the  current  and  previous  values,  and  the 
process  is  continued  until  the  temperature  into  the  partial 
node  exactly  equals  the  saturation  temperature.  The  values 

for  Q,  T  X.  ,  T-  .  and  U  .  are  stored  to  be  used  later  to 

'  g  out  f  in  oi 

develop  weighted  averages  with  the  boiling  portion  parameters. 

The  remainder  of  the  node  is  solved  in  the  boiling 
portion  of  the  model. 


E.  BOILING  NODAL  ANALYSIS 

Nodes  in  the  boiling  section  are  solved  by  the  same 
iterative  scheme  developed  in  the  superheating  section.  The 


correlations  utilized  to  determine  the  solution  are  differ¬ 
ent  since  flow  in  the  boiling  regime  will  involve  two  phases 
on  the  fluid-side.  In  two  phase  flow  the  inside  heat  trans¬ 
fer  coefficient  will  change  with  quality,  x,  and  the  heat 
flux,  q".  The  correlation  selected  for  this  model  is  the 
seune  as  that  used  by  Combs  [Ref.  3]  which  was  recommended 
by  Tong  [Ref.  7]  for  both  nucleate  boiling  and  forced  convec 
tion.  The  amalysis  presented  by  Combs  is  presented  here. 

The  equation  for  the  two-phase  flow  heat  transfer  coeffi 
cient  was  given  by  Schrock  and  Grossman  in  Tong  [Ref.  7]  as 
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The  constants  are  given  by  Wright  [Ref.  7]:  B  =  6.70  x  lo  , 
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A  »  3.5  X 10  ,  and  n  =  0.66.  The  heat  transfer  coefficient, 

assxuning  a  totally  liquid  flow  is 


d^  G(l-x)  0.8  C_,  0.1 


0.023  3^[ 


and  the  Martinelli  parameter  is 
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0.9  p  0.5  y  0.1 
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The  limits  of  the  data  for  this  correlation  a_3; 


quality:  0.05 — 0.57 


heat  flux:  6.0x10^  —  1.45x10^  BTU/hr-ft^, 


In  this  model,  heat  flux  in  the  boiling  section  will  be  at 
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levels  considerably  below  6  x lo  BTU/hr-ft  and  the  full 
range  of  quality  must  be  considered.  Therefore,  in  order  to 
predict  the  performance  of  the  Schrock  and  Grossman  correla¬ 
tion  at  lower  heat  flux  and  quality  edjove  0.57,  a  plot  was 
made  (Fig.  7),  of  the  ratio  where  h^^  is 

the  heat  transfer  coefficient  from  the  Dittus-Boelter  corre¬ 


lation  for  water  at  saturated  liquid  conditions.  The  corr'^- 
lation  of  Dengler  and  Addoms,  as  given  in  Tong  [Ref.  7], 
indicates  that  low  heat  flux  and  high  water- vapor  mixture 
velocity  favor  the  forced  convection  mechanism.  After  a 
review  of  the  information  available  in  Tong,  it  is  not 
entirely  clear  whether  the  forced  convection  or  the  nucleate 
boiling  mechanism  is  dominant  for  the  relatively  low  flow 
velocities  and  low  heat  flux  of  this  model.  In  any  case, 
the  Schrock  and  Grossman  correlation  applied  at  low  heat 
flux  represents  a  type  of  compromise  between  the  entirely 
forced  convection  assumptions  of  Dengler  and  Addoms  and  the 
"mixed"  assumption  of  Schrock  and  Grossman. 

A  simplification  assumed  in  the  model  is  that  the 
quality  in  the  node  is  the  average  quality  for  the  entire 
length  of  the  node.  Also,  to  expedite  the  numerical  proce¬ 
dure,  it  is  assumed  that  the  node  effectiveness  is  0.25, 
initially.  The  effectiveness  rarely  attains  a  value  higher 
than  this  for  this  type  heat  exchanger,  so  the  assumption  is 
valid.  The  value  of  the  effectiveness  is  adjusted  to  the 
"true"  value  in  the  iterative  process.  In  a  process  involving 
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a  change  of  phase  the  fluid  is  isothermal.  In  actuality, 
there  is  a  slight  temperature  change  through  the  boiling 
regime  but  for  purposes  of  analysis  this  change  is  neglected. 
For  this  reason  the  fluid  bulk  temperature  in  any  node  in 
the  boiling  portion  is  equal  to  the  fluid  saturation  tempera¬ 
ture.  Therefore,  to  determine  the  amount  of  energy  cJosorbed 
by  the  fluid  in  a  node,  the  change  in  fluid  enthalpy  between 
the  inlet  and  outlet  is  calculated. 

In  the  first  node  in  the  boiler,  which  is  the  node  imme¬ 
diately  before  the  superheater  or  a  partial  node  if  the  node 
was  split,  the  gas  bulk  temperature  is  again  set  equal  to 
the  gas  inlet  temperature.  The  gas-side  heat  transfer  coeffi 
cient  is  calculated  using  the  same  procedure  as  delineated 
in  the  superheater  nodal  analysis.  Next,  a  Q  for  the  node 
is  calculated  from 


Q  =  e  m  C  (T  .  -  T,  .  ) 

g  Pg  y  in  f  in 


where  the  effectiveness,  e,  is  initially  assumed  to  be  0.25 
and  T^  is  equal  to  the  fluid  saturation  temperature.  A 
fluid  inlat  enthalpy  is  derived  from 


'^f  in  ^f  satv  "  ^ 


where  h^  satv  *  enthalpy  of  the  fluid  at  saturated  vapor. 


The  fluid  bulk  enthalpy  is 
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=  (h-  .  +  h_ 
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The  average  quality  for  the  node  is  calculated  from 

X  =  ~  satl 

^f  satv  "  ^f  satl 

Tidiere  h^  satl  enthalpy  of  the  fluid  at  a  saturated 

liquid  condition.  The  Martinelli  parameter  is  determined 
from 
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where  P^/P^  the  liquid/vapor  densities  at  the  saturation 
temperature  and  are  the  liquid/ vapor  viscosities  at 

the  saturation  temperature.  Now  the  heat  transfer  coefficient 
for  totally  liquid  flow  is  calculated 


h^  =  0.023  gi 
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From  this  the  heat  tremsfer  coefficient  is  calculated  for 
two-phase  flow 
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where  h,^  *  h  -  „  . .  -  h_  ,  and  the  constants  A,  B  and  n 

rg  x  satv  t  satl 

have  the  aforementioned  values  and  G  is  the  mass  flow  rate 


per  unit  area.  Finally,  the  node  overall  heat  transfer 
coefficient  is  calculated  from 
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Whenever  two  phase  flow  is  involved,  the  heat  capacity  of 
the  fluid  is  assumed  infinite.  For  this  reason,  within 
the  boiling  regime,  will  always  equal  the  heat  capacity 

of  the  hot  side,  so. 
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and  the  Ntamber  of  Transfer  Units,  NTU,  is 
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The  node  effectiveness  is 


e  =  1  -  exp ( -NTU)  . 


Now,  a  revised  Q  for  the  node  is  calculated  from 
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From  this  Q  a  revised  gas  outlet  temperature  is  calculated 
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With  this  gas  temperature  out,  a  new  film  temperature  is 
calculated  in  the  same  manner  as  for  the  superheating  section. 

The  gas  outlet  temperature  is  stored  and  the  film  temp¬ 
erature  is  used  to  repeat  the  process  for  the  node.  At  the 
end  of  the  second  time  through,  the  gas  outlet  temperature 
derived  is  compared  with  the  result  from  the  first.  If 
equality  exists,  a  solution  for  the  node  has  been  found  and 
the  model  proceeds  to  the  next  node.  If  inequality  exists, 
the  result  of  the  current  iteration  is  stored  for  comparison 
with  the  results  of  the  following  iteration. 

When  a  final  solution  is  attained,  the  enthalpy  of  the 
fluid  at  the  inlet  to  the  node  is  checked  against  the 
enthalpy  of  saturated  water  at  the  heat  exchanger  pressure. 

If  the  inlet  enthalpy  is  greater  than  or  equal  to  the  enthalpy 
of  saturated  water,  the  model  proceeds  to  the  next  node  and 
begins  the  solution  process  again.  If  the  inlet  enthalpy 
is  less  than  the  enthalpy  of  saturated  water  then  the  model 
proceeds  to  a  node  splitting  routine  similar  to  that  used 
when  transitting  from  the  superheating  regime  to  the  boiling 
regime . 

The  node  splitting  is  again  accomplished  by  halving  the 
mass  flow  rate  on  the  gas-side  and  halving  the  physical 
length,  L,  of  the  node.  The  remaining  geometric  parameters 
are  recomputed  and  the  solution  process  is  reinitiated.  The 
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mass  and  length  adjustment  procedure  is  continued  until  the 

inlet  enthalpy  is  equal  to  the  enthalpy  of  saturated  water. 

When  this  is  achieved,  the  values  of  T  ,  U  .  and  Q  are 
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stored  to  be  recombined  as  weighted  averages  with  the  heating 
portion  of  the  node  to  give  the  overall  nodal  values. 

The  remainder  of  the  node  is  solved  in  the  heating  portion 
of  the  model. 

F.  HEATING  NODAL  ANALYSIS 

The  procedure  for  analyzing  a  node  in  the  heating  portion 
of  the  model  is  identical  to  that  used  in  the  superheating 
portion.  Each  node  is  analyzed  in  an  iterative  process,  and 
it  makes  no  difference  if  the  first  node  was  a  split  node. 

The  first  node  in  the  heating  section  will  have  a  gas  inlet 
temperature  equal  to  the  gas  outlet  temperature  of  the  node 
directly  below.  The  fluid  outlet  temperature  will  equal  the 
fluid  saturation  temperature  for  the  specified  pressure. 

An  outline  of  the  calculations  is  as  follows: 

1.  Calculate  the  gas-side  heat  transfer  coefficient 
using  the  node  gas  inlet  temperature  as  the  gas  bulk 
temperature . 

2.  Calculate  the  fluid-side  heat  transfer  coefficient 
initially  using  the  fluid  saturation  temperature  as 
the  fluid  bulk  temperature. 

3.  Calculate  the  overall  heat  transfer  coefficient  using 
the  effectiveness-NTU  method. 


9 .  store  the  result  of  the  first  pass  through  the  model 
and  iterate  until  procedure  converges  to  solution, 
utilizing  gas  outlet  temperature  for  comparison  at 
each  pass . 

All  properties  of  the  gas  and  fluid  are  evaluated  at  the 
respective  bulk  temperatures  for  the  node.  The  process  is 
repeated  node  by  node  until  the  node  at  the  fluid  inlet  to 
the  heat  exchanger  is  completed.  The  final  node  is  either 
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node  (1,1)  for  an  odd  number  of  passes,  or  node  (1,10)  for 
an  even  number  of  passes  in  the  heat  exchanger. 

G.  HEAT  EXCHANGER  SOLUTION  PROCEDURE 

When  the  procedure  of  the  heating  section  has  arrived  at 
a  solution  for  the  final  node  of  the  heat  exchanger,  the 
fluid  inlet  temperature  for  the  node  is  compared  to  the  user 
specified  feedwater  inlet  temperature.  If  equality  exists, 
the  model  proceeds  to  an  output  routine  where  the  matrices 
of  water/steam  temperatures,  exhaust  gas  temperatures, 
overall  heat  tremsfer  coefficients  and  individual  heat 
transfer  per  node  are  printed.  Additionally,  the  fluid-side 
mass  flow  rate,  minimum  pinch  point  temperature  difference, 
total  Q  for  the  heat  exchanger  and  isentropic  steam  turbine 
power  are  calculated  and  printed. 

The  pinch  point  temperature  difference  is  defined  in 
this  model  as  the  minimum  difference  in  the  heat  exchanger 
between  the  gas  inlet  temperature  and  the  fluid  inlet  temp¬ 
erature  to  any  node.  Consideration  of  the  effect  of  a  pinch 
point  below  approximately  25  deg.  F  is  cause  for  concern 
only  if  the  low  pinch  point  is  coupled  with  a  low  gas  temp¬ 
erature.  This  combination  is  non-prohibitive  to  hot  side 
corrosion.  If  the  final  node  fluid  inlet  temperature  is  not 
equal  to  the  user  specified  feedwater  inlet  temperature,  then 
the  fluid  mass  flow  rate  that  was  "guessed"  at  the  beginning 
is  adjusted  and  the  entire  model  is  run  again.  If  the  solution 


temperature  is  too  high  the  "guessed"  fluid  mass  flow  rate 
is  lowered.  With  less  mass  going  through  the  heat  exchanger 
the  flow  velocity  is  decreased.  A  lower  flow  velocity  through 
the  heat  exchanger  enables  more  energy  from  the  gas-side  to 
be  transferred  to  the  fluid-side.  A  larger  Q  will  increase 
the  temperature  differential  between  the  fluid  inlet  and 
outlet  which  will  yield  a  lower  solution  fluid  inlet  temp¬ 
erature  the  second  time  through  the  model.  The  process  is 
opposite,  i.e.,  mass  is  added,  if  the  fluid  inlet  temperature 
to  the  first  node  is  found  to  be  too  low. 

The  model  will  continue  to  adjust  the  mass  after  each 
iteration  in  smaller  amounts  until  an  exact  match  is  made  to 
the  user  specified  feedwater  inlet  temperature.  When  equality 
is  achieved,  the  model  proceeds  to  the  output  routine  where 
the  adaove  mentioned  parameters  are  printed. 

Included,  as  Appendix  A,  are  sample  outputs  of  the  model 
for  the  uniform  distribution,  distribution  1  and  distribu¬ 
tion  5  for  all  four  of  the  selected  gas  mass  flow  rates. 

The  run  numbers  at  the  top  of  the  first  page  of  each  output 
were  utilized  for  accovinting  purposes.  The  uniform  distri¬ 
bution  was  used  for  r\ins  1,  10,  19  and  28;  distribution  1  was 
used  for  rvins  2,  11,  20  and  29;  and  distribution  5  was  used 
for  runs  6,  15,  24  and  33.  The  first  page  of  the  output  shows 
the  inter-active  questions  asked  by  the  model  and  the  user 
responses.  The  following  pages  show  the  matrices  developed 
by  the  model  in  the  solution  process .  Comparison  of  the 
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matrices  for  the  uniform  and  non-uniform  distributions 
demonstrates  the  serious  impact  of  the  non- uniformity  on  the 
performance.  The  non-uniform  distribution  profiles  can  be 
seen  as  projected  profiles  on  the  results  of  each  of  the 
matrices.  The  final  page  of  the  output  is  a  listing  of 
additional  parameters  which  are  calculated  by  the  model  includ¬ 
ing  the  fluid  mass  flow  rate,  Q  based  on  the  individual 
matrices,  location  of  the  beginning  of  boiling  and  super¬ 
heating  within  the  WHRU  and  the  isentropic  horsepower. 


III.  DISCUSSION  OF  RESULTS 

A.  INITIAL  TESTING 

The  foregoing  model  was  coded  in  FORTRAN  IV  and  applied 
in  a  computer  simulation  on. an  IBM-3033  computer. 

To  properly  test  the  model  and  develop  an  understanding 
of  the  operation  of  the  WHRU  several  runs  were  made  utilizing 
a  uniform  gas  distribution  with  a  total  mass  flow  rate  and 
temperature  corresponding  to  a  gas  turbine  horsepower  rating 
for  approximately  a  20  kt.  speed.  The  conditions  presupposed 
usage  of  a  General  Electric  LM2500  Gas  Turbine  Engine,  and 
gas  mass  and  temperature  data  was  taken  from  [Ref.  8] .  Initially 
the  model  was  run  with  attention  to  pinch  point  temperatures. 

The  ntimber  of  passes  was  varied  to  find  the  largest  heat  ex¬ 
changer  which  could  be  used  and  still  maintain  a  pinch  point 
in  excess  of  25  degrees  fahrenheit.  Under  conditions  of  non- 
uniform  flow  distributions  the  model  definition  of  pinch  point 
is  misleading  since  small  gas  mass  flow  rates  through  a 
particular  gas  path  may  yield  a  nearly  zero  pinch  point  when 
the  overall  performance  of  the  heat  exchanger  is  not  affected. 
Inspection  of  the  model  output  for  the  non-uniform  distribu¬ 
tions  will  demonstrate  this  result.  The  preliminary  runs 
allowed  selection  of  a  23  pass  heat  exchanger  with  a  29.4  deg. 

F  pinch  point.  Therefore,  the  23  pass  heat  exchanger  became 
the  "design."  The  performance  data  tabulated  for  the  LM2500 
engine  [Ref.  8]  was  utilized  to  select  four  gas  mass  flow 
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rates  and  inlet  temperatures  which  correspond  closely  to 
data  points  on  a  cubic  power  curve.  The  four  combinations 
of  mass  and  temperature  chosen  were: 


”g  ' 

121.8 

Ibm/s , 

T  . 
g  in 

*  844®  F, 

SHP  =  17500 

,  NPT  =  3000  rpm 
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99.1 

Ibm/s , 

T  . 
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=  755“  F, 

SHP  =  10000, 

NPT  =  2400  rpm 

*9  ' 

77.4 

Ibm/s , 

T 
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=  702“  F, 

SHP  =  5000, 

NPT  =  1800  rpm 

m  = 

g 

66.7 

Ibm/s , 

T 

9  in 

=  694“  F, 

SHP  =  3000, 

NPT  =  1200  rpm 

Each  of  these  was  run  in  the  model  as  uniform  and  various 
non-uniform  distributions. 

To  develop  the  non-uniform  distributions,  a  normalized 
flow  distribution  table  was  constructed  and  is  presented  as 
Table  I.  This  table  shows  the  fraction  of  the  total  mass 
flow  rate  which  is  passing  through  each  individual  gas  path 
for  each  of  nine  different  flow  distributions.  The  flow 
distributions  were  arbitrarily  derived  in  an  attempt  to 
develop  a  spectrum  of  distributions  that  may  enter  the  heat 
exchanger.  A  graphical  display  of  each  of  the  flow  distribu 
tions  is  shown  in  Figures  8A-8I.  The  total  mass  flow  rate 
used  was  multiplied  by  each  entry  in  the  normalized  flow 
distribution  table  to  arrive  at  the  mass  flow  rate  through 
each  of  the  ten  gas  paths  for  any  of  the  nine  distributions. 
This  procedure  allowed  a  comparison  of  the  effects  of  non- 
uniform  distributions  on  the  performance  of  the  heat  ex¬ 
changer.  All  nine  distributions  were  run  for  each  of  the 
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four  chosen  mass  and  temperature  combinations.  The  actual 
mass  flow  distributions  used  are  shown  in  Tables  II-A — II-D. 
The  combination  of  nine  distributions  and  four  different 
total  gas  flow  rates  yielded  36  runs.  For  the  runs,  the 
remaining  parameters  were: 


a) 

feedwater  inlet  temperature 

= 

200 

deg. 

F 

b) 

desired  steam  temperature 

= 

650 

deg. 

F 

c) 

steam  operating  pressure 

= 

400 

psia 

d) 

minimum  allowed  gas  outlet  temperature 

= 

300 

deg. 

F 

e) 

node  length 

1.0 

ft . 

f) 

niimber  of  tubes  per  row 

= 

32 

g) 

geometric  scaling  factor 

= 

1.0 

The  model  calculates  the  isentropic  horsepower  that  would 
be  developed  by  an  ideal  steam  turbine.  This  power  can  also 
be  viewed  as  the  maximum  power  available  from  the’ developed 
superheated  steam  exhausting  to  a  fixed  condenser  pressure. 

The  method  of  calculation  is  to  determine  the  enthalpy  of  the 
steam  at  the  outlet  pressure  and  temperature  and  then  sub¬ 
tract  the  enthalpy  of  the  steam  after  isentropic  expansion 
to  a  condenser  pressure  of  2  psia.  The  difference  in  enthalpy 
was  multiplied  by  the  fluid  mass  flow  rate,  with  the  appro¬ 
priate  conversions,  to  arrive  at  a  horsepower  value.  The 
horsepower  calculated  for  each  run  is  tabulated  in  Table  III. 
The  t£U3le  lists  the  results  of  each  run  of  the  model  for  all 
36  runs  for  each  of  5  different  steam  outlet  temperatures. 
Table  IV  is  a  similar  table  showing  the  solution  values  of 
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the  fluid  mass  flow  rates  derived  in  each  run.  The  two 
tables  demonstrate  the  variation  of  fluid  mass  flow  rate  and 
power  which  is  attained  by  varying  the  desired  steam  outlet 
temperature;  and  the  effect  of  the  non-uniform  distributions 
on  performance. 

B.  EFFECT  OF  NON-UNIFORM  DISTRIBUTIONS 

The  data  from  column  1  of  Table  III  was  plotted  and  is 
shown  in  Figures  9A-9E  as  power  vs.  gas  flow  rate.  Figure 
9A  shows  the  comparison  of  the  uniform  distribution  with 
distributions  1  and  2.  Figure  9B  shows  the  comparison  of 
the  uniform  distribution  with  distributions  3  and  4,  etc. 

From  the  graphs  it  can  be  seen  that  the  distribution  inflicting 
the  most  serious  impact  on  the  heat  exchanger  performance  is 
distribution  1  (see  Figure  8-B) .  Distribution  1  is  highly 
center  concentrated  and  utilizes  the  area  available  for  heat 
treuisfer  the  least  efficiently.  For  comparison,  distribution 
3  is  a  nearly  uniform  distribution  and  it  inf licts  a  relatively 
minor  impact  on  the  performance.  Distributions  2  and  4  are 
ranked  next  in  order  of  impact.  These  are  also  distributions 
where  the  mass  is  centrally  concentrated,  but  to  a  lesser 
degree  than  distribution  1.  Distributions  5  and  6  are  mirror 
images  of  one  another  and  also  have  a  zero  mass  flow  rate 
through  two  of  the  gas  paths.  A  zero  mass  flow  rate  indicates 
that  the  flow  through  that  particular  pass  has  no  velocity. 

In  a  real  situation  there  will  be  a  transient  heat  transfer 
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from  the  still  mass  of  gas  surrounding  the  tubes  by  natural 
convection  until  the  gas  temperature  reaches  the  tube  wall 
or  fluid  temperature.  In  any  case,  the  amount  of  heat 
transferred  in  such  a  node  is  negligible  compared  to  the 
nodes  experiencing  a  gas  flow  rate.  For  the  model,  a  zero 
gas  flow  rate  in  a  node  will  be  treated  as  having  zero  heat 
transfer,  i.e.,  a  "dead  node,"  therefore,  the  Q  for  the  node 
was  set  to  zero,  the  fluid  inlet  temperature  and  fluid  outlet 
temperature  were  set  equal  and  the  gas  outlet  temperature 
was  set  equal  to  the  gas  inlet  temperature.  In  reality,  the 
gas  temperature  in  a  "dead  node"  will  attain  the  temperature 
of  the  fluid  in  the  node  but  for  pinch  point  calculation  pur¬ 
poses  the  original  gas  inlet  temperatiare  was  maintained. 

Also,  since  the  model  assumes  no  cross-mixing  of  the  gas 
between  gas  paths  as  the  distribution  passes  through  the 
heat  exchanger,  this  assumption  ignores  the  potential  tendency 
of  the  gas  mass  distribution  to  become  uniform  during  passage. 
Distributions  5  and  6  are  skewed  to  the  left  and  right, 
respectively.  The  combination  of  the  null  gas  flow  paths  and 
the  skewedness  of  these  distributions  caused  the  second  most 
detrimental  effect  on  the  heat  exchanger.  It  should  be  noted, 
however,  that  whether  the  distribution  was  skewed  to  the  right 
or  left  side  had  relatively  no  effect;  that  is,  the  results 
for  distribution  5  are  nearly  identical  to  those  for  distri¬ 
bution  6.  The  same  result  was  true  for  distributions  7  and 
8,  which  are  also  mirror  images  of  each  other  with  left  and 
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right  skewedness.  The  effect  of  distributions  7  and  8  were 
very  nearly  the  same  as  those  for  distribution  4.  Figure  9E 
is  a  ccxnposite  which  demonstrates  the  effect  of  all  the  dis¬ 
tributions.  The  ranking  as  to  effect  yields: 

1.  the  uniform  distribution  gives  the  best  performance 

2.  distr.  3,  nearly  uniform 

3.  distr.  4,  least  center  concentrated,  nearly  same  as 
distr. 's  7  and  8 

4.  distr.  2,  intermediately  center  concentrated 

5.  distr. 's  5  and  6,  skewed  with  null  flow  paths 

6.  distr.  1,  highly  center  concentrated 

The  ranking  can  be  viewed  as  a  measure  of  the  efficiency 
of  the  distribution  in  utilizing  the  area  available  in  the 
heat  exchanger  for  heat  transfer.  As  the  amount  of  distor¬ 
tion  of  the  flow  increases,  the  quantity  of  energy  trans¬ 
ferred  decreases. 

Figures  lOA-lOE  are  graphical  representations  of  the 
solution  fluid  mass  flow  rates  vs.  the  gas  mass  flow  rates 
for  each  of  the  distributions;  exactly  as  was  done  to  develop 
Figures  9A-9E.  As  was  to  be  expected,  the  results  for  Figures 
lOA-lOE  are  identical  to  those  for  power  vs.  gas  flow  rate, 
since  power  is  directly  proportional  to  the  fluid  mass  flow 
rate. 

After  these  preliminary  results,  it  was  decided  to  con¬ 
duct  the  remaining  runs  utilizing  only  the  uniform  distribu¬ 
tions  and  the  worst  of  the  non-uniform  distributions;  i.e., 


distribution  1.  It  was  assimed  that  results  for  the  remain¬ 
ing  distributions  would  fall  between  the  values  of  these  two 

C.  EFFECT  OF  HEAT  EXCHANGER  SIZE 

Next,  an  analysis  was  done  to  determine  the  effect  on 
output  power  and  fluid  mass  flow  rate  by  varying  the  nvimber 
of  heat  exchanger  passes.  Using  the  four  chosen  gas  mass 
flow  rates,  runs  were  conducted  for  heat  exchangers  of  20, 
21,  23  and  25  passes.  All  runs  were  for  650  deg.  F  outlet 
steam  at  400  psia  operating  pressure.  Therefore,  the  runs 
for  the  23  pass  heat  exchanger  were  the  "design"  runs  and 
the  runs  for  20,  21  and  25  passes  were  the  "off-design"  r\ans 
provided  for  comparison.  The  results  of  these  runs  are 
tabulated  in  Table  V  and  graphically  presented  in  Figures 
llA-llH  as  power  and  fluid  mass  flow  rate,  respectively,  vs. 
the  number  of  heat  exchanger  passes.  Again,  because  of  the 
relationship  between  power  and  fluid  mass  flow  rate,  the 
graphical  results  are  similar  in  contour.  The  results  show 
that  as  the  size  of  the  heat  exchanger  is  increased,  so  will 
the  power  out.  To  provide  the  increased  power  the  mass  flow 
rate  of  the  fluid  must  increase  since  for  the  different 
sized  heat  exchangers  the  inlet  and  outlet  enthalpies  of  the 
fluid  is  not  changing.  The  limit  to  the  increase  in  size  is 
the  available  energy  in  the  gas.  Increasing  the  heat  ex¬ 
changer  size  to  the  point  where  sufficient  energy  is  trans¬ 
ferred  from  the  gas  to  the  fluid  to  attain  the  minimum 


acceptable  exhaust  gas  temperature  is  the  maximiim  size 
possible. 

The  curves  display  an  interesting  result  concerning  the 
differences  between  the  uniform  distribution  and  the  non- 
uniform  distribution.  As  the  size  of  the  heat  exchanger 
increases,  the  two  curves  are  converging  with  the  difference 
in  power  (or  fluid  mass  flow  rate)  being  less  for  the  larger 
heat  exchanger  than  for  the  smaller.  From  the  data  of  Table 
V,  it  can  be  seen  that,  for  the  121.8  Ibm/s  gas  flow  rate, 
the  power  from  distribution  1  is  89%  of  the  power  attained 
from  the  uniform  distribution  for  a  20  pass  heat  exchanger. 
This  increases  to  approximately  91%  for  the  25  pass  heat 
exchanger.  The  more  rapid  increase  in  power  for  the  non- 
\iniform  distribution  is  caused  by  the  mass  concentration  in 
the  center  gas  paths.  For ‘a  particular  size  heat  exchanger, 
the  paths  with  a  small  gas  mass  flow  will  attain  the  temper¬ 
ature  of  the  fluid,  usually  in  the  boiling  region.  From  that 
point  on,  these  paths  contribute  little,  or  not  at  all,  to 
the  total  heat  transfer.  The  high  concentration  of  mass 
in  the  center  paths,  however,  contribute  greatly  to  the  total 
heat  transfer  because  of  the  high  velocity  of  the  flow  and 
the  high  heat  transfer  coefficients.  These  paths  of  high 
concentration  will  contribute  until  the  gas  exits  the  heat 
exchanger.  By  adding  more  passes  to  the  heat  exchanger 
it  is  as  if  the  "dead  areas"  in  the  low  mass  flow  paths  were 
placed  in  the  path  of  the  high  mass  flows.  The  effect  is  a 


significant  increase  in  the  heat  transfer.  The  high  heat 
treuisfer  coefficients,  mass  velocity  and  temperature  differ¬ 
entials  all  contribute  to  a  greater  percentage  increase  per 
pass  than  is  achieved  by  the  uniform  distribution.  As  addi¬ 
tional  passes  are  added,  the  performance  of  the  non-uniform 
distribution  improves,  but  it  will  never  attain  the  level  of 
the  uniform  distribution.  For  two  heat  exchangers  of  equal 
size,  the  uniform  distribution  gives  the  most  efficient  use 
of  the  available  area  for  heat  transfer. 

D.  EFFECT  OF  STEAM  OUTLET  TEMPERATURE 

Each  distribution,  uniform  and  distributions  1-8,  was 
run  at  each  of  the  four  gas  flow  rates  for  five  specified 
steam  outlet  temperatures.  The  power  and  fluid  mass  flow 
calculated  by  the  model  are  presented  in  Tables  III  and  IV. 

The  results  are  also  presented  graphically  in  Figures  12-A 
and  12-B  for  the  uniform  distribution  and  distribution  1 
at  each  of  the  four  gas  flow  rates. 

Both  the  tabular  and  graphical  data  shows  a  nearly  uni¬ 
form  increase  in  fluid  mass  flow  rate  with  the  linear  de¬ 
crease  in  steam  outlet  temperature.  The  tabular  data  for 
power  demonstrates  a  non-linear  result,  showing  initially  a 
decrease  in  power  with  a  decrease  in  steam  outlet  temperature, 
and  then  increasing  to  attain  a  power  higher  than  that  for 
a  650  deg.  F  steam  outlet  temperature. 

The  non-linearity  of  the  power  is  caused  by  the  non¬ 
linearity  of  the  steam  enthalpy  as  the  steaun  outlet  temperature 
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decreases.  The  combination  of  the  change  in  enthalpy  and 
the  change  in  the  fluid  mass  flow  rate  results  in  a  mislead^ 
ing  impression  of  attaining  a  higher  power  at  a  lower  steam 
temperature.  The  total  power  does  increase  but  the  specific 
power  decreases.  The  increase  in  fluid  mass  flow  between  a 
650  deg.  F  steam  outlet  temperature  and  a  550  deg.  F  steam 
outlet  temperature  is  approximately  7%.  The  increase  in 
power  at  the  lower  steam  outlet  temperature  is  less  than  0.2%. 
The  penalty  to  utilize  a  lower  steam  temperature  emd  main¬ 
tain  an  equivalent  power  output  is  a  water/steam  system  which 
must  be  able  to  circulate  7%  more  fluid  in  the  same  time 
period.  This  result  may  have  serious  impact  in  the  design 
when  condenser  size,  circulating  pump  size,  valves,  piping 
and  the  steam  turbine  are  taken  into  consideration.  The 
steam  turbine  must  also  be  considered  with  respect  to  the 
lower  steeun  temperature  having  a  lower  quality  through  the 
latter  stages  of  the  expansion.  The  effects  of  overly  "wet" 
steam  can  be  disastrous  to  many  turbines.  The  end  result 
is  that  the  sacrifice  in  power  at  the  higher  steam  tempera¬ 
ture  is  justified  by  the  decrease  in  fluid  flow  required 
and  higher  end  point  quality. 

E.  EFFECT  OF  HEAT  EXCHANGER  PRESSURE 

The  uniform  distribution  and  distribution  1  were  run  for 
each  of  the  four  gas  flow  rates  with  three  different  speci¬ 
fied  heat  exchanger  operating  pressures.  The  pressures  were 
the  "design"  pressure  of  400  psia  and  then  500  psia  and  600 


psia.  The  tabulated  results  of  these  runs  are  presented  in 
Table  VI,  and  graphically  presented  in  Figures  13-A  and  13-B 
for  power  and  fluid  mass  flow  rate,  respectively,  vs.  heat 
exchanger  pressure.  Again,  because  of  the  direct  coLzsla- 
tion  between  fluid  mass  flow  rate  and  power,  the  contours  of 
the  graphical  results  are  similar. 

The  curves  show  a  decrease  in  power  and  fluid  mass  flow 
rate  as  heat  exchanger  pressure  is  increased.  As  the  pres¬ 
sure  is  increased  by  50%  above  the  design  level,  there  is  a 
4%  drop  in  che  system  fluid  mass  flow  rate  and  a  7.5%  loss 
in  power.  Clearly,  the  benefits  derived  from  the  decrease 
in  system  fluid  mass  flow  rate  does  not  justify  the  sacri¬ 
fice  in  system  power.  With  respect  to  the  analysis  of  the 
previous  section,  a  compromise  must  be  reached  to  achieve 
the  maximum  power  from  the  system  within  the  confines  of  a 
maximum  acceptable  fluid  mass  flow  rate. 

The  graphical  results  also  show  the  effect  of  the  non- 
uniform  distribution  for  each  of  the  four  gas  flow  rates. 

As  the  gas  flow  rate  decreases,  the  penalty  for  non-uniformity 
has  a  lesser  effect.  The  percentage  loss  for  the  four  flow 
rates  are  approximately  9,  8,  7  and  6  percent,  respectively, 
from  highest  flow  rate  to  lowest.  The  loss  is  approaching 
zero  as  the  flow  rate  approaches  zero.  The  loss  in  power  is 
directly  related  to  the  decrease  in  fluid  mass  flow  rate 
determined  in  the  solution;  and  the  decrease  in  fluid  flow 
rate  is  a  result  of  the  ineffectiveness  of  the  non- uniformity 
on  the  operation  of  the  heat  exchanger. 


F.  EFFECT  OF  GAS  FLOW  RATE  ON  HEAT  EXCHANGER  PERFORMANCE 

The  data  from  the  matrix  of  water/steam  temperatures  for 
the  four  gas  flow  rates  as  uniform  distributions  was  plotted 
as  a  function  of  tube  length.  The  curves  are  shown  in 
Figures  14-A — 14-D.  Figxire  14-E  is  a  plot  of  the  water/ 
steam  temperatures  as  a  function  of  tube  length  for  distri¬ 
bution  1  with  a  gas  mass  flow  rate  of  121.8  Ibm/s.  Figure 
14-F  is  a  superposition  of  Figure  14-E  on  Figure  14-A. 

Analysis  of  Figures  14-A  through  14-D  shows  the  effect 
of  the  decrease  in  gas  flow  rate  on  the  rate  of  heat  trans¬ 
fer  as  the  fluid  passes  through  the  heat  exchanger.  The 
higher  the  gas  flow,  the  steeper  is  the  slope  of  the  curve 
in  the  superheat  region.  A  steeper  slope  is  indicative  of 
a  shorter  tube  length  required  to  attain  the  desired  super¬ 
heat  temperature.  The  effect,  as  the  gas  flow  rate  is 
decreased,  is  that  more  of  the  heat  exchanger  is  dedicated 
to  superheating.  The  length  of  the  tube  required  to  raise 
the  water  from  a  saturated  liquid  to  a  saturated  vapor 
condition  is  nearly  the  same  for  all  gas  flow  rates  with 
only  a  "left  shift"  on  the  plots  caused  by  the  superheat 
conditions.  The  slope  of  the  curves  in  the  heating  region 
is  the  opposite  of  that  of  the  superheat.  As  the  gas  flow 
rate  decreases  the  slope  of  the  curve  in  the  heating  region 
increases. 

Since,  for  all  conditions,  the  feedwater  inlet  tempera¬ 
ture  and  steam  outlet  temperature  are  the  same,  all  tubes 


are  continuous  and  of  the  scune  length,  the  difference  in  the 
curves  is  attributable  to  the  rate  of  heat  transfer  and  the 
total  amount  of  energy  transferred.  The  rate  equation  is 


Q  =  C  (T  .  -  T 
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The  change  in  gas  mass  flow  rate  will  obviously  decrease  the 
total  amount  of  energy  transferred  since  less  gas  in  means 
less  energy  into  the  heat  exchanger.  If  the  curves  were 
plotted  in  a  third  dimension  with  Q,  the  smallest  Q  would 
occur  with  the  smallest  gas  flow  rate.  The  forced  equality 
with  the  rate  equation  means  that,  as  the  gas  flow  rate 
decreases,  the  rate  of  heat  transfer  must  decrease.  The 
figures  show  the  total  effect  on  the  profile  operation  of 
the  heat  exchanger  is  of  little  consequence.  The  similarity 
of  the  curves  is  an  indication  of  the  "physical  operation" 
of  the  heat  exchanger  vinder  varying  conditions.  For  the  test 
r\ins,  a  decrease  in  gas  flow  rate  was  coupled  with  a  decrease 
in  gas  inlet  temperature.  The  two  parameters  changing  and 
the  similarity  of  the  curves  emphasizes  the  physical  constancy 
of  the  design. 


Figure  14-E  is  the  plot  for  distribution  1.  The  "ripple" 
effect  in  the  heating  and  superheating  regions  is  a  result 
of  the  non- uniformity  of  the  gas  mass  distribution,  but  the 
superposition  shown  in  Figure  14-F  demonstrates  the  oblivious¬ 
ness  of  the  physical  properties  of  the  heat  exchanger  to  a 
variation  of  input. 

G.  EFFECT  OF  NON-UNIFORM  MASS  DISTRIBUTION  ON  EXHAUST  GAS 

TEMPERATURE 

For  all  runs  of  the  model  the  temperature  distribution 
of  the  exhaust  gas  entering  the  heat  exchanger  was  assumed 
uniform.  To  detemine  the  effects  of  a  non-uniform  gas  mass 
distribution  on  the  exhaust  gas  temperature  profile,  a  plot 
was  made  of  the  gas  temperatures  across  several  of  the  heat 
exchanger  passes.  The  matrix  of  exhaust  gas  temperatures 
was  used  and  data  was  plotted  for  both  a  uniform  distribu¬ 
tion  and  distribution  1.  The  gas  mass  flow  rate  was  121.8 
Ibm/s.  Figures  15-A  and  15-B  show  the  results. 

Figure  15-A  demonstrates  that  with  a  uniform  mass  and 
temperature  distribution  at  the  inlet  to  the  heat  exchanger, 
the  distribution  is  virtually  undistorted  by  passage  through 
the  WHRU.  This  result  must  be  considered  ideal  with  respect 
to  the  assumption  of  no  cross-mixing  of  the  gas  during  passage. 
There  is  a  smooth  transition  of  energy  from  inlet  to  outlet 
which  accounts  for  the  best  performance  of  the  heat  exchanger 
under  this  condition.  The  effect  of  the  non- uniformity  is 
clearly  seen  in  Figure  15-B.  The  temperature  distribution  is 


62 


rapidly  distorted  to  become  almost  a  replica  of  the  gas  flow 
distribution.  This  is  a  result  of  the  disparity  of  heat 
transferred  and  heat  transfer  rates  from  node  to  node  in 
the  WHRU.  A  slight  smoothing  effect  is  seen  in  the  boiling 
portion,  but  distortion  is  again  accentuated  in  the  heating 
section.  In  a  real  situation  there  will  be  some  non¬ 
uniformity  of  temperature  as  it  enters  the  heat  exchanger. 

The  analysis  shows  that  the  condition  will  be  worsened  to 
the  degree  of  the  mass  non-uniformity  with  the  subsequent 
loss  of  performance  in  the  WHRU. 

The  distortion  of  the  exhaust  gas  temperature  distribu¬ 
tion  forces  consideration  of  pinch  point.  When  the  tempera¬ 
ture  differential  between  the  exhaust  gas  and  the  tube 
becomes  less  than  approximately  25  deg.  F  coupled  with  a  gas 
temperature  less  than  approximately  300  deg.  FV  the  WHRU  may 
experience  localized  hot-side  corrosion  due  to  sulphur  products 
in  the  exhaust  gas.  It  should  be  noted  that  this  condition 
will  exist  only  near  the  top  of  the  heat  exchanger  near  the 
exhaust  gas  exit.  If  the  degree  of  flow  distortion  is 
severe  and  maintained  for  a  long  period,  the  effects  can  be 
disastrous . 
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IV.  CONCLUSIONS 


The  waste  heat  recovery  unit  is  a  practical  method  of 
recouping  thermal  energy  from  gas  turbine  engine  exhaust 
gases  for  fuel  economy  improvement.  In  the  design  and 
installation  of  a  RACEK  system  several  factors  must  be  taken 
into  consideration. 

1.  Mass  Flow  Distribution.  The  performance  of  a  waste 
heat  recovery  unit  is  adversely  affected  by  a  non-uniform  gas 
mass  distribution;  and  the  degradation  is  directly  propor¬ 
tional  to  the  degree  of  non-uniformity.  To  ensure  acceptable 
results,  the  system  should  have  flow  distribution  control 
devices  installed  between  the  gas  turbine  and  the  WHRU. 

2.  Operating  Pressure.  The  operating  pressure  of  the 
system  affects  the  power  output  inversely.  A  higher  pressure 
system  requires  less  fluid  to  deliver  the  same  power  but  the 
percent  reduction  in  fluid  required  ^oes  not  justify  the 
percent  loss  in  power.  Since  the  system  is  not  a  main  pro¬ 
pulsion  system  requiring  a  large  volume  of  fluid,  the  model 
suggests  using  a  low  pressure  system. 

3.  Steam  Outlet  Temperature.  Reducing  the  outlet  temp¬ 
erature  required  forces  an  increase  in  fluid  mass  flow  rate 
to  maintain  the  same  power  out.  High  temperature  steam,  how¬ 
ever,  has  a  greater  wear  effect  on  system  components  such  as 
valves,  seals  and  piping.  The  model  suggests  a  steam  tempera- 
tvire  as  high  as  practical  should  be  used,  within  material 
constraints . 


4.  Heat  Exchanger  Size.  Power  and  fluid  mass  flow  rate 
are  directly  related  to  heat  exchanger  size.  The  size  of 
the  WHRU  should  be  as  large  as  the  available  space  and 
weight  constraints  allow,  or  as  large  as  necessary  to  pro¬ 
duce  the  power  required.  Consideration  must  be  given  to  the 
system  as  a  whole  since,  as  the  heat  exchanger  becomes  larger, 
so  does  the  piping,  valves,  pumps  and  condenser  necessary  to 
operate . 

5.  Speed.  Power  developed  is  a  function  of  the  energy 
available  in  the  exhaust  gas  flow.  At  higher  gas  turbine 
speeds  more  horsepower  is  availcdsle  from  the  WHRU  to  assist 
the  main  propulsion  system.  Utilizing  the  latest  NAVSEA 
speed  profile,  the  heat  exchanger  should  be  sized  in  accordance 
with  the  ship  speed  most  commonly  used. 

All  data  suggests  a  WHRU  as  large  as  physically  possi¬ 
ble,  operated  at  a  low  steam  pressure,  high  steam  temperature, 
and  with  flow  distribution  control  devices  installed  provides 
the  best  perf ormcuice . 


RECOMMENDAT IONS  FOR  FURTHER  RESEARCH 


A.  WHRU  MODEL  IMPROVEMENT  AND  EXPANSION 

The  control  imposed  on  the  WHRU  model  in  this  study  was 
adjustment  of  the  steam  flow  rate  to  maintain  a  user  speci¬ 
fied  steam  outlet  temperature.  This  forces  the  steam  system 
to  respond  to  a  variation  of  gas  inlet  conditions  by  "throt¬ 
tling**  the  feedwater  system  to  maintain  the  required  steam 
temperature . 

This  is  only  one  of  several  types  of  control  systems 
which  may  be  used  in  a  RACER  system  installation.  Two 
possible  alternative  systems  cure:  one  which  maintains  a 
constant  steam  flow  rate  while  allowing  the  steam  tempera¬ 
ture  to  fluctuate,  and  one  in  which  pressure  is  allowed  to 
"float"  with  steam  demand.  In  this  second  type  system,  the 
flow  resistance  of  the  system  is  used,  and  as  the  fluid  mass 
flow  rate  decreases  (steam  demand  decreases) ,  the  pressure 
required  decreases.  Knowing  the  inter-relationship  will 
yield  a  "floating"  as  well.  Investigation  into  the 

performance  of  a  WHRU  with  these  types  of  control  could  lead 
to  an  improved  output.  The  investigation  of  all  types  of 
control  should  consider  the  effects  on  the  entire  system 
because,  the  type  of  control  may  dictate  the  type  steam 
turbine  which  must  be  used. 

To  expand  the  understanding  of  the  effects  of  non-uniform 
gas  distributions  on  the  WHRU  performance,  the  current 


model  could  be  enlarged  to  allow  a  variation  of  the  input  in 
an  additional  dimension.  A  two  dimensional  non-uniform  input 
will  yield  a  three  dimensional  picture  of  the  temperature 
distributions  as  the  gas  passes  through  the  heat  exchanger. 

An  analysis  of  this  type  would  help  design  the  flow  distri¬ 
bution  control  devices  to  eliminate  severe  trouble  spots 
contributing  to  hot  side  corrosion  and  to  maximize  the  WHRU 
performance . 

In  addition  to  the  expansion  to  a  three  dimensional 
analysis,  the  current  model  could  be  improved  by  including 
the  following  items. 

1.  Fluid-Side  Pressure  Drop 

The  water/s team- side  pressure  drop  may  be  a  signifi¬ 
cant  percentage  of  the  system  operating  pressure.  Calculation 
of  the  pressure  drop  and  its  effect  on  system  performance 
should  be  included. 

2.  Additional  Geometries 

Additional  fin- tube  configurations  should  be  con¬ 
sidered  for  possible  improvement  of  the  heat  transfer  charac¬ 
teristics  and  the  weight  and  space  requirements  of  the  WHRU. 

3.  Fouling 

Consideration  of  inside  and  outside  heat  transfer 
surface  fouling  should  be  given  when  investigating  possible 


fin-tvibe  configurations. 

4 .  Alternate  Working  Fluids 


a  minimxiia  cost  in  system  maintainability,  while  maintaining 
safety  to  operating  personnel. 

5.  Non-Uniform  Temperature  Distributions 

The  current  model  could  be  utilized  to  determine  the 
effects  of  non-uniform  temperature  distributions.  The  gas 
is  at  a  uniform  temperature  as  it  exits  the  gas  turbine,  but 
it  probably  is  not  uniform  at  the  inlet  to  the  WHRU.  Inves¬ 
tigation  into  the  effects  of  non-uniformity  of  both  gas  mass 
and  temperature  should  be  accomplished. 

The  model  should  be  incorporated  into  a  complete 
RACER  system  cuialysis,  including  models  of  the  steam  turbine, 
condenser  and  circulating  pump.  Variations  of  gas  distribu¬ 
tion  input  can  then  be  analyzed  to  determine  total  system 
effect.  Understanding  the  total  system  effect  and  the  inter¬ 
relationship  of  the  varicLbles  could  then  be  used  to  optimize 
the  system. 


TABLE  II-A 


GAS  FLOW  DISTRIBUTIONS  (ihg  =  121.8  LBM/S) 


CAS  MASS  PLOW  MATE  ■  121 <8  L9M/S 

CAS  -tEMMEMAIUllE  AT  heat ‘CXCHANGEM  INLET  ■  8«A  DEC.  P 
CASE  STS  tWEPS  I  ■  irSOO  NAT  m  3009 


12.18  I  12.18  !  12.18  I  12.18 


12.18  12.18 


S.MT  I  8.28  1  13. S*  !  38.9*  I  30.94  I  13.64  !  8.28  |  S.9T  )  2.87 


3.78  I  S.38  I  IO.«7  I  18.20  t  25.09  |  25.09  I  l4.20  I  10. *7  S.38  3.78 


10.98  I  I  2.10  I  13.48  I  14.82  114.82  )  13.40  |  12.18  |  10.98  |  9.74 


S.05  I  7.92  I  12. 79  17.05  1  20.10  1  20.  10  1  17. OS  |  12.79  7.02  |  3. OS 


27.28  I  23.29  19.95  t  15.96  i  11.69  7.80  |  3. 65  |  0.80  |  0.00 


7.80  I  11.69  IS. 96  19. 8S  |  23.29  |  27.28  |12.I8 


iS.tO  j  10.47  8.69  I  8.11 


5.24 


IDBIE&BQnQDI 


I  S.  47 


TABLE  II-B 


GAS  FLOW  DISTRIBUTIONS  (itig  =  99.1  LBM/S) 

CAS  MASS  Plow  rate  >  99.1  lsm/s 

CAS  TEMOCRAfUOS  AT  MEAT  EXChANSET  1 N^ET  s  7S5  DEC.  P 
CASS  561  I  REP  81  5mP  «  10000  nPT  b  2400 


NPt  B  2400 


6.74  I  I . to  2S.I 7  2S. I  7  I  11.10 


e.S2  I  13. 18  20.41  I  20.41  13.18 


10.90  11.90  11.90  I  10.90 


8.44  I  10.41  I  19.37  18.95  16.35 


9.91  I  22.20  I  19*03 


1  2.59  19*23  12.29 


16.19  I  12.99 


9.5l  I  12.99  I  16.15 


8.50  7.23 


19.23  12.59 
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TABLE  II-C 


GAS  FLOW  DISTRIBUTIONS  (m^  ^  11 A  LBM/S) 


«AS  mtk%%  FLO*  KATe  ■  rr.A  lbm/s 

6AS  TexvcBArtfae  at  hcat  eacmangc*  tNkCT  «  T02  oec*  r 
CASE  SS3  fSEFS  I  ShB  •  SOOtt  M»T  ■  ISOO 


a.o«  I  o«oo 


9.83  lS.02 


6AS  VATM 
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19. 6S  19.44 


4.4S  10.29  I  IS. 94  I  15.94  |  10.29 
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5.03 


7.74  47.3a  14.04 


2.  32 
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7. A3  4.94 


7.43  10. IS  12. 4t 


4.4S 


5.45  I  4.45 


o 
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10 

■89 

7.  7* 

7.74 

S.24 

3.00 

I.3I 

4.45 

3.41 

2.41 

■SI 

4.97 

0.  1  3 

S.03 

1.94 

2.32 

0.00 

0.00 

14.04 

1  7.3* 

7.74 

4. 93 

3.33 

1  12.15 

tS.02 

9.83 

TABIE  II-D 

GAS  FLOW  DISTRIBUTIONS  (m  =66.7  IWS) 

9 

GAS  MASS  7L3W  RATE  «  67.7  LBM/S 

GAS  TEMMERATiiftE  AT  hEA T  EXCHANGER  INLET  «  69A  OEG.  * 

CASE  945  I  REF  8)  SHP  m  3000  M»T  «  1200 


BIST. 

NO. 


CAS  PATH 

3  4 


4.77 


3.32  I  A. 40 


2.90 


4.10 


4.40  I  7.11 


S.  I  7  I  3. 00 


2.03 


4.77 

1  4.77 

4.77 
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1  7.20 
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i 
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7.44 
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0.  13 

■BB 

11.17 
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0.00 

4.  SO 

4.34 
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5.02 
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4.50  1 

4.9* 

5.02 

0.00 

0.40 

2.91 


11.03 

4.50 


I  0.43 


13.1* 


4.77 

2.91 


TABLE  III 


K-.' 


ISENTROPIC  POWER  FOR  VARIOUS  GAS  FLOW  DISTRIBUTIONS, 
GAS  FLOW  RATES  AND  STEAM  OUTLET  TEMPERATURJ S 


RUN 

NO. 

OISTR. 

NO. 

STEAM  TEMPERATURE 

650  625  600  575  550 

GAS 

MASS 

121.8 

LBM/S 

1 

U 

6944 

6943 

6933 

6946 

6956 

2 

1 

6286 

6315 

6295 

6317 

6348 

3 

2 

6550 

6570 

6550 

6580 

6598 

4 

3 

6924 

6925 

6910 

6927 

6939 

5 

4 

6692 

6708 

6691 
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6  730 

6 

5 

639  7 

6419 

6415 

6420 

6434 

7 

6 

6390 

6382 

6407 

6431 

6464 

8 

7 

6754 

6742 
~  6733 

6722 

6727 

6750 

9 

8 
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6744 

6763 

6762 
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MASS 
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10 

U 
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11 

1 
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4183  1 

12 

2 
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4323 

13 

3 

44/5 
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14 

4 

4  3S6 
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15 

5 
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16 

6 
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4194  ;  4187 

4220 

4245  1 

17 

7 
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18 

8 
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4417 
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MASS 

77.4 
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19 

U 
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2957 
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20 

1 
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21 

2 
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2872 

2632 

22 

3 

2943 
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2  962 

2  970 

23 

4 

2882 
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24 

S 
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2  833 

2837 

25 

6 
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2811 
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26 

7 
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2914 
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27 

B 
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2933 

GAS 

MASS 

67.7 
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28 

U 
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2522 

2  530 
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29 

1 

2365 
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23  93 

24U 
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30 

2 
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2471 

31 

3 
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2527 

2533 

32 

4 

246  5 
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2480 

2490 

2498 

33 

5 

2377 

240  5 

2418 

2435 

2437 

34 

6 

2388 

2406 

2416 

2425 

2442 

35 

7 

2476 

2481 

2492 

2496 

2499 

36 

8 

246  9 

2487 

2482 

2496 
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TABLE  IV 

FLUID  MASS  FLOW  RATE  FOR  VARIOUS  GAS  FLOW  RATES,  GAS  FLOW 
DISTRIBUTIONS  AND  STEAM  OUTLET  TEMPERATURES 


GAS 

MASS 

121.6 

LBM/S 


S 

SS 
.1 
LBH/S 


RUN 

OISTR. 

STEAM  temperature 

NO. 

NO. 

650 

625 

600 

575 

550 

B 

U 

12.33 

12.54 

12.75 

12.97 

13.20 

2 

1 

11.17 

11.41 

11.58 

1 1 . 79 

12.04 

3 

2 

11.63 

11.87 

12.05 

12.23 

12.52 

12.30 


1 1.6  9 


11.36 


12.00 


12.51 


12.12 


11.59 


11.53 


12.16 


6.11 


7.46 


2 


12.31 
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12 .53 
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12.40 

!  12.63 

8.25 

8  .40 

7.64 
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13.17 


12.77 


12.21 


12.83 


S.56 


7.94 


8.20 


8.54 


8 


3 


8.05 


8 


8.38 


5 


5.33 


5 


5 


5 


5.29 

1  5.38 

5.27 

i  5.39 

5.45 

5.54 

5.45 

5.56 

4.72  i 

4.81 

TABLE  V.  POWER  AND  FLUID  FLOW  RATE 

FOR  VARIOUS  SIZE  HEAT  EXCHANGERS 


NUMBER  OF  HEAT  EXCHANGER  PASSES  UNITS 


TABLE  VI.  POWER  AND  FLUID  FLOW  RATE  FOR  VARIOUS 
STEAH  OPERATING  PRESSURES 


DISTR. 

OPERATING  PRESSURE 

400  500  600 

UNITS 

■■ 

U 
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HP 

1 

6286 

KM 

HP 

BIH 

121.8 

U 

12.33 
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1 
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u 
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HP 
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1 
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HP 
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1 
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HP 
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u 

MBM 

2339 

Bi 

GAS 

1 

2  208 

MASS 

BH 

1 

Bi 
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u 

A. 23 

4.01 

LBM/S 

1 

BB 

3  .99 

3.78 
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LAYOUT  OF  MARINE  GAS  TURBINE  WASTE-HEAT  STEAM  GENERATOR:  ONCE-THROUGH, 
COUNTER-CROSS  FLOW  HEAT  EXCHANGER 


Figure  1.  Waste  Heat  Recovery  Unit  Layout 


NR 


Nuinber  of  Heat  Exchanger  Passes 


General  Node  (m,n} 


Figure  2.  Heat  Exchanger  Nodal  Arrangement 
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Figure  5.  Tube  Fin  Configuration 
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Figure  6.  Waste  Heat  Recovery  Unit  T\ibe  Layout 
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Figure  7.  vs.  Quality 


GRS  FLOW  DISTRIBUTION 
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Figure  8-B.  Flow  Distribution  #1 


FLOW  DISTRIBUTION 


Figure  8-D.  Flow  Distribution  #3 
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S  FLOW  DISTRIBUTION 


Figure  8-F.  Flow  Distribution  #5 


S  FLOW  DISTRIBUTION 


Figure  8-H.  Flow  Distribution  #7 


Figure  8-1.  Flow  Distribution  #8 


POWER  VS.  GflS  MASS  FLOW  RATE 


Figure  9-A.  Power  vs.  Gas  Plow  Rate  (Distr 
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GRS  MASS  FLOW  RRTF  (LBM/S) 


RRTE  (LBM/S) 


POWER  VS.  GPS  MRSS  FLOW  RRTE 


Figure  9-D.  Power  vs.  Gas  Flow  Rate  (Dlstr 
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Fluid  Flow  Rate  vs.  Gas  Flow  Rate  (Distr 


FLUID  FLOW  RRTE  VS.  GRS  FLOW  RRTE 


RRTE  CLBM/S) 


POWER 


Figure  11-B.  Power  vs.  No.  of  Heat  Exchanger  Passes  (lUg  =  99.1  Ibm/s) 


•yr**.** 


Figure  11-D.  Power  vs.  No.  of  Heat  Exchanger  Passes  (m  =66.7  Ibm/s) 


Figure  11-E.  Fluid  Flow  Rate  vs.  No.  of  Heat  Exchanger  Passes 
(ra  =  121.8  Ibm/s) 


flSS  FLOW  RRTE  VS.  NO.  OF  PRSSES 


Fluid  Flow  Rate  vs.  No.  of  '  ’  •»nger  Passes 

(m  =99.1  Ibm/s) 


Figure  11-G.  Fluid  Flow  Rate  vs.  No.  of  Heat  Exchanger  Passes 
(m  =77.4  Ibm/s) 


STEAM  OUTLET  TEMP  (DEG  F) 


FLUID  FLOW  RBTE  VS.  STERM  TEMP 


Figure  12-B.  Fluid  Mass  Flow  Rate  vs.  Steam  Outlet  Temperature 
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HEAT  EXCHANGER  PRESSURE  (PSIA) 


Figure  13-B.  Fluid  Mass  Flow  Rate  vs.  Hea 
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Figure  14-A.  Fluid  Temperature  vs.  Tube  Length  (m  =  121.8  Ibm/s) 


Figure  14-B.  Fluid  Temperature  vs.  Tube  Length  (m  =  99.1  Ibm/ 
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Figure  14-F.  Fluid  Temperature  vs.  Tube  Length 
(Distr.'s  U,  1;  m  =  121.8  Ibra/s) 


TEMPERRTURE  DISTRIBUTION 


Figure  15-A.  Gas  Temperature  Distribution  (Uniform,  121.8  Ibm/s) 


APPENDIX  A 


Sample  runs  of  the  computer  output  are  provided  for  a 
comparison  of  the  effects  of  the  non-uniform  gas  mass 
distributions  on  the  temperatures  within  the  heat  exchanger 
and  the  effect  on  overall  performance. 

Runs  1/  2  and  6  utilized  a  total  gas  flow  rate  of 
121.8  Ibm/s  with  the  uniform  distribution,  distribution  1  and 
distribution  5,  respectively.  Runs  10,  11  and  15  utilized 
the  same  distributions  but  with  a  total  gas  flow  rate  of 
99.1  Ibm/s.  Runs  19,  20  and  24;  and  runs  28,  29  and  33 
also  utilized  the  same  distributions  but  with  mass  flow 
rates  of  77.4  Ibm/s  and  66.7  Ibm/s,  respectively. 
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